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ABSTRACT 

We present a study of the prevalence of optical and radio nuclear activity with respect 
to the environment and interactions in a sample of SDSS galaxies. The aim is to deter- 
mine the independent effects of distinct aspects of source environment on the triggering 
of different types of nuclear activity. We defined a local density parameter and a tidal 
forces estimator and used a cluster richness estimator from the literature to trace dif- 
ferent aspects of environment and interaction. The possible correlations between the 
environmental parameters were removed using a principal component analysis. By far 
the strongest trend found for the active galactic nuclei (AGN) fractions, of all AGN 
types, is with galaxy mass. We therefore applied a stratified statistical method that 
takes into account the effect of possible confounding factors like the galaxy mass. We 
found that (at fixed mass) the prevalence of optical AGN is a factor 2-3 lower in the 
densest environments, but increases by a factor of ~ 2 in the presence of strong one- 
on-one interactions. These effects are even more pronounced for star-forming nuclei. 
The importance of galaxy interactions decreases from star-forming nuclei to Seyferts 
to LINERS to passive galaxies, in accordance with previous suggestions of an evolu- 
tionary time-sequence. The fraction of radio AGN increases very strongly (by nearly 
an order of magnitude) towards denser environments, and is also enhanced by galaxy 
interactions. Overall, the results agree with a scenario in which the mechanisms of 
accretion into the black hole are determined by the presence and nature of a supply 
of gas, which in turn is controlled by the local density of galaxies and their interac- 
tions. A plentiful cold gas supply is required to trigger star-formation, optical AGN 
and radiatively-efficient radio AGN. This is less common in the cold-gas-poor environ- 
ments of groups and clusters, but is enhanced by one-on-one interactions which result 
in the flow of gas into nuclear regions; these two factors compete against each other. In 
the denser environments where cold gas is rare, cooling hot gas can supply the nucleus 
at a sufficient rate to fuel low-luminosity radiatively-inefficient radio AGN. However, 
the increased prevalence of these AGN in interacting galaxies suggests that this is not 
the only mechanism by which radiatively-inefficient AGN can be triggered. 

Key words: galaxies: evolution - galaxies: interactions - galaxies: active - radio 
continuum: galaxies - catalogues 



1 INTRODUCTION 

Active galactic nuclei (AGN) are closely related to galaxy 
formation and evolution. The black holes that power AGN 
are found in all massive galaxies and their masses are 
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tightly correlated with both the masses and the veloc- 
ity dispersions of the stellar bulges (e.g. Marconi fc Huntl 



120031 : iFerrarese fc MerrittfeOOOl ; iGebhardt et alj|2000l ). AGN 
may play an important role in the feedback mechanisms 
th at control the growth of massive galaxies (see review 
by ICattaneo et al.l |2009i , and references therein). Interac- 
tion between galaxies can drive gas i nto the central region 
of the galaxy and trigger the AGN (Shlosman et al.l [l990l ; 
iBarnes fc Hernquisth99ll ; lHaan et al.ll2009l ; lLiu et alfeoill ). 
However, a hi gh fraction of AGN may by fueled by secu lar 
processes (e.g. ISilverman et al.|[20lTI : (Sabater et al]|2012i ). 
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The effect of the large-scale environment and of galaxy 
interactions on the triggering of an AGN has previously 
been studied by many authors, but contradictory results 
have been found. Quas ars are thought to be associated 



with g alaxy interact i ons ( Sanders et al.|[l988l ; lUrrutia et al.l 



120081 : iLetawe et al ] 12010) or a higher number of com- 
panions dHutchings fc Campbelllll983l ; iBahcall et al] 1 19971 ; 



ISerber^al] M- s lv^allt^dirf ound a direct ^ 



AGN (Petrosianl 1982 


; Koulouridis et al.l 20061; Alonso et al. 


2007; iRogers et al. 




20091: Ellison et al.l 2011; Liu et al. 


20121: iHwang et alJ 


2012). On the other hand, other stud- 



ies d id not find this relation 
200 ll: 



IIBushousel 119861; iSchmittl 



Ellison et al] 120081; liTet all 120081 ; iDarg et al.l l20ld ; 



Slavcheva-Mihova fc Mihovl l201ll ). In some cases, a lower 



prevalence of AGN is found in environments with higher 
density of galaxies dCarter et all l200ll ; iMiller et ail 120031 ; 
IKauffmann et al.1 l20oJ )7 or the AGN prevalence is found 
to depend upon distance to the centre o f the nearest 
group or clu s ter l|yon der Linden et al.ll2010l) . In contrast, 
lAlonso et al.l (|2012T ) found an increase of the fraction of pow- 
erful AGN toward denser environments in a sample face- 
on spiral galaxies. Some studies find signa t ures of mergers 



in ac t ive galaxies (IComerford et al.l 120091 ; ISmirnova et al 
|2010|; IVillar-Martm et al.l l201ll , I20121 ). but. iGrogin eta! 
1 20051 ') found no connection between recent gal axy merg- 



ers an d A GN. Making use o f radio observations, iKuo et al.l 
( 2008) and lTang et al. |2008l) measured the disruption of the 
atomic gas distribution and kinematics to trace interactions 
and found a clear relation with the presence of an AGN. 

With respect to X-ray sel ected AGN, a hig her incidence 
of signs of disruption in AGN (Koss et al. 2010) and a higher 



fraction of AGN in pairs dSilverman et al " 120111 ) and in the 
centre of clusters l|Ruderman fc Ebelind 120051 ) have been 
found. However, a relation between in t eraction and AGN 
was n ot found by Pierce et al. d2007l). Georgakakis et al.l 
(|2009h or lGabor et alJ (|2009l ). IWaskett et al.1 (|2005l ) found 
that the environments of AGN are indisti nguishable from 
those of normal, inactive galaxies. Finally, ISilverman et al.l 
l|2009l ) found no difference between the fraction of AGN in 
groups and in field galaxies. 

If radio selected AGN are considered, iBest et all 
( 2005a]) found a tendency for radio-loud AGN t o be l ocated 
in richer environments. iReviglio fc Helfandl (2006) con- 
firmed this trend for radio AGN, which is in contrast to that 
widely found for optical AGN. Brightest group and cluster 
galaxies are more likely to host a radio-lou d AGN than othe r 
galaxies with similar s tellar masses (e.g. IBest et al.l [20071) . 
iDomingue et al.l l|2005l ) found an excess of radio AGN in 
pairs of g alaxies while no rad i o AG N were f ound in isolated 
galaxi es (jSabater et al.1 120081 . l20ld ) . While iHeckman et al.l 
(1986) found many ra dio AGN t o be "strongly peculiar" 
in optical morphology, iTal et all l|2009l ) found no correla- 
tion between radio AGN and tidal distortions. Neverthe- 
less, there is recent evidence of the triggering of powerful 
radio galaxies by in teractions |Ramos Almeida et~aT1l2012l ; 
iBessiere et ai1l2012l ). 

Part of the confusion from the results discussed above 
may involve the different definitions of "environment" and 
"interaction" used in the previous studies. This may be: a) 
a higher local density of neighbours, b) the membership of 
a cluster or group versus field or filaments, c) relative loca- 



tion within a group or cluster, or, d) one-on-one interactions. 
These definitions can lead to contradictory results since they 
trace different aspects of interaction and environment. It is 
known that many interacting g alaxy pairs c an be found in 
low density regions (even void; IVerlev et al.l 2007; Argudo- 
Fernandez et al. in prep.), consequently, the concept "field" 
is not equivalent to "isolated" . We might expect the amount 
and physical conditions of the gas in galaxies that is used to 
fuel the AGN to b e related to large-scale environment (e.g. 
ISiiacki et al1l2007l ). On the other hand, interactions at the 
one-on-one level are expected to set the gravitational dis- 
turbances that could trigger the AGN. Therefore, it is fun- 
damental to separate those different aspects when studying 
their effect on nuclear activity. 

A further complication is the type, and luminosity, 
of the AGN selected in different studies. Although the 
emission of an AG N spans a wide range of the electro- 
magnetic spectrum (|Elvis et al.l 1199-3 ). the different meth- 
ods and wavelengths used to select AGN will also lead to 
the selection of different kinds of AGN. All of the selec- 
tion methods are affected by selection effects (iHd 2008). 
for example, X-rays have been used to detect AGN tha t 
did not present optical emission lines ( Mar tini et al.l [20021 ). 
Selection based on optical emission lines can discern be- 
tween Seyfert and low-ioniza tion nuclear emission-line re- 
gion (LINER; lHeckmanlll980l ) types, however, the final clas- 
sification may depend on the exact processing method used 
and the classification criteria a p plied (e.g. iKewlev et al.l 
l2006l ; IConstantin fc Vogelevl2006l ). rBest et al.l (|2005al ) found 
that the probability of galaxy harbouring a radio AGN 
is independent of its optical classification and conclude 
that low-luminosity radio AGN and emission-line AGN are 
powered by different physical processes. Two main feed- 
ing mechanisms can be found in AGN: a) "standard" high- 
excitation radiatively-efficient mode, and, b) low-excit ation 
radiatively-inefficient mode fsee lBest fc Heckman|[2012l . and 
references therein). Radiatively-efficient AGN include lumi- 
nous radio AGN (also called High Excitation Radio Galaxies 
or HERG) and the optical or X-ray AGN while radiatively- 
inefficient ones are observed as low-lumin osity radio AGN 
(Low Excitation Radio Galaxies or LERG; lHardcastle et al.l 
120071). These types of feeding m echanisms can evolve in time 
( Hlavacek-Larrondo et al.ll2012l ) and may be fundamental to 
exp lain the decline of star formation in large elliptical galax- 
ies dCattaneo et al.ll2009l ). 

It is also important to consider that many properties 
of galaxies depend on the environment, with clusters typ- 
ically hosting mor e massive and more early-type galaxies 
than the field (e.g. Blanton fc Mou stakas 2009; Den g et"ail 
120091 ; IPark fc Choill2009l ; iDeng et all l201ll ). If samples are 
not well-matched (as has been the case for many previous 
studies) then these differences can give rise to false depen- 
dencies. In particular it is known that galaxy mass is one 
of the most (if not the most) im portant driving factors fo r 
the preva lence of: a) optica l (e.g. IKauffmann et aill2003bl), 
b) radio (IBest et alJ 20o"5al ). and c) X-ray ( Silverman et al.1 
120091 ; iTasse et al.ll201ll ) selected AGN. The p revalence of 
AGN also depends upon gal a xy morphology dMoles et al.l 
Il995l ; ISchawinski et al.1 |201(]| ; ISabater et all |2012| ) These 
properties are also closely related to the environment via 
the density-morphology and density-luminosity relations 
(|Dresslerlll980l ; IKauffmann et al.ll2004l : iBlanton et ai1l2005al ; 
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iDeng et al.ll2009l; IPark fc Choill2009l ; iDeng et al.ll201ll ). In- 
deed, Park fc Choil l 2009^ suggest that galaxy properties are 
mainly driven by mass and morphology, and environment 
plays only a secondary role. Hence, these confounding fac- 
tors should be carefully taken into account in order to obtain 
unbiased results. 

We aim to study independently the effects of the large- 
scale environment and the smaller-scale interactions on the 
triggering of both radio and optically selected AGN. In 
Sect. [21 the sample and the data used are presented. We 
define a set of parameters to quantify the environment and 
the interactions that affect our galaxies, and these are de- 
scribed in Sect. [3] The results of the statistical studies that 
estimate the relation between the environment and inter- 
action with the prevalence of different types of active nu- 
clei are presented in Sect. [4] Finally, these results are dis- 
cussed in Sect. [5] and the final conclusions are presented 
Sect. [6] Throughout the paper, the following cosmologi- 
cal parameters are assumed: f2 m = 0.3, J1a = 0.7 and 
H = TOkms^Mpc" 1 . 



2 THE SAMPLE AND THE DATA 
2.1 The sample 

The sample was based on the seventh data release (DR7 
lAbazaiian et al. | 20091) of the Sloan Digital Sky Survey 



(SDSS; York et alj 2000T ). The base sample is composed 



of galaxies from the main s pectroscopic sample (magnitudes 
between 14.5 < r < 17.77; IStrauss et"ai]|2002h with a red- 
shift between 0.03 and 0.1. Galaxies with a redshift below 
0.03 were discarded, a) to limit the search for companions to 
a reasonable area of the sky around the target galaxies, b) 
to lessen the possible effect of the incompleteness of galaxies 
with high brightness and, c) to avoid the possible errors in 
the measurements of the photometry due to the large size 
of the galaxies. Galaxies with a redshift above 0.1 were also 
discarded, a) to allow a good signal-to-noise for the emission 
lines used for the optical activity classification and for the 
radio-continuum emission used in the radio activity classifi- 
cation, and, b) to maximize the number of companions with 
similar luminosities that are covered by the spectroscopic 
survey. The evolutionary time between redshift 0.03 and 0.1 
is w 9 x 10 s years We take no account of any evolution of 
the nuclear activity prevalence during this period, but this is 
supposed to be small and in any case any evolutionary effects 
are minimized by the use of the stratified statistical study. 
The sample is composed of galaxies from the table Galaxy 
which have an assigned spectroscopic object specObjID with 
a spectroscopic redshift (field bsz) within the limits. 

A small percentage of luminous galaxies are not in- 
cluded in the main spectroscopic survey at lower redshifts 
(those with r < 14.5 which at z = 0.03 corresponds to 
L r ~ 1O 1O ' 2 L ). At z = 0.1 the magnitude limit corre- 
sponds to a luminosity limit in r-band of ~ 1O 9,4 L0. For 
most of our analyses, stellar mass (see next Section) limits 
of 10 ^ log(M/M0) 12 were applied to reduce the in- 
completeness. In any case, in our study we compute relative 
fractions of galaxies instead of absolute numbers. These frac- 
tions should not be biased by the small number of missing 
galaxies, especially since we observe no significant trends of 



any of the AGN fractions with redshift across our samples. 
The use of relative fractions, the cut in masses, and the use 
of stratified statistical methods that take into account the 
mass of galaxies should mitigate the possible effect of the 
incompleteness. 

The environmental parameters that we aimed to obtain 
could be affected if the target galaxy is close to the border 
of the SDSS footprint or near to a bright star. The poten- 
tial loss of a significant fraction of companion galaxies could 
then cause an underestimation of the parameters. Hence, 
we discarded from our study galaxies that could be affected 
by these problems. To do this we used the New York Uni- 
versity (NYU) Value-Added Galaxy Catalog (|Blanton et al.l 
2005b), which provides additional data to the SDSS sample. 
The parameter NEXPECT found in the simpleden table gives 
an estimation of the number of tracer galaxies that can be 
found within a given fixed radius of a target galaxy. These 
tracer galaxies are homogeneously and randomly distributed 
but only within the coverage region of the SDSS spectro- 
scopic survey. Hence, if the number is lower than a given 
value, it indicates that the region around the target galaxy 
is not completely covered by the survey. We discarded galax- 
ies with NEXPECT < 0.45. In some cases, there were no avail- 
able data from the NYU catalogue and we took the value of 
NEXPECT from the closest galaxy with data within a maxi- 
mum of 30 arcmi n radius. A vi sual inspection (using the soft- 
ware TOPCAT; lTavloiil2005l ) showed that this process was 
reliable for the selection of galaxies in well-sampled SDSS 
areas. Four galaxies were manually added and three regions 
close to the borders (126 galaxies in total) were manually 
discarded. Two additional regions containing 1594 galaxies 
close to the borders were nagged later after finding system- 
atic offsets for those galaxies in one of the environmental 
parameter that we computed (the tidal estimator; Sect. 13. 2] ). 

The final sample was composed of 267977 galaxies. Re- 
liable estimations of the cluster richness using a friends of 
friends algorithm (as explained in Sect. 13. 3[) could only be 
obtained for galaxies brighter than M r — 5 log(h) = —20. 
These galaxies constitute the reduced sample which was used 
when the cluster richness parameter was needed for the anal- 
ysis. The total number of galaxies in each sample is shown 
in Table Q] 

2.2 AGN classification 

Our aim was to check the dependence of the different nuclear 
activity types on the environment and interaction. For this 
study, we considered optically and radio selected nuclear 
activity types. 

Data based on optical spectra were drawn from the 
Max Plank Institute for Astrophysics and Johns Hopkins 
Uni versity (MPA-JHU) adde d value spectroscopic catalogue 
(cf. iBrinchmann et al.ll2004h. Information about total stel- 
lar mass l|Kauffmann et al.ll2003ah and the classification of 
the optical nuclear activity based on the accurate measure- 
ment of emission-line fluxes can be found in this catalogue. 
The classification of optical activity is based on the typica l 
diagnostic diagrams (BP T diagrams; iBaldwin et"al] Il98ll ; 
iKauffmann et alj l2003br i. The following optically-selected 
nuclear activity types are distinguished: a) star forming nu- 
clei (SFN), b) transition objects (TO; which may contain a 
mix of star-formation and an AGN) c) Seyfert, d) LINER, 
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Table 1. Number of galaxies of each type. 



Type_^ Whole Reduced Whole Reduced 

sample sample sample sample 

M lim£ M lim. 

Total 267977 107176 201425 106503 



Passive 99180 47333 90174 47311 

Optical AGN 57955 34266 55864 34255 

LINER 16288 12315 16144 12313 

Seyfert 7170 3691 6763 3689 

TO 34497 18260 32957 18253 

SFN 110842 25577 55387 24937 



Radio AGN 1137 1015 1129 1015 

LERG 1087 979 1082 979 

HERG 47 34 44 34 



a Meaning of the different types: Passive - galaxies with 
^[Oiii] < 10 6,5 Lq; Optical AGN - galaxies with L[om] ^ 
10 6 ' 5 Lq classified as Low Ionization Nuclear Emission Re- 
gions (LINER), Seyferts or Transition Objects (TO); SFN - 
Star Forming Nuclei galaxies with Lrom] > 10 6 - 5 L ; Radio 
AGN - galaxies with Li.4GHz > 10 3 W m -2 Hz -1 classi- 
fied as radio AGN including High Excitation Radio Galaxies 
(HERG) and Low Excitation Radio Galaxies (LERG). 
b Mass limited to 10 s: log(M/M Q ) ^ 12 



2.3 Potential companions 

For our tidal forces estimator (see Sect. I3.2|l . when inves- 
tigating potential companions for our targets we need to 
consider not only the SDSS spectroscopic galaxies but also 
galaxies with similar luminosities that were not covered by 
the spectroscopic survey (e.g. due to fibre allocation restric- 
tions) and galaxies with lower luminosities (below the spec- 
troscopic magnitude limit) that could be companions of our 
targets. We considered the objects in the SDSS Galaxy cata- 
logue with: a) r magnitude below 22.0, and b) spectroscopic 
redshift below 0.11, or photometric redshift minus three 
times its error below 0.10. With this magnitude limit we in- 
clude potential companions at least « 50 less luminous than 
the least luminous target galaxy at a given redshift. The 
spectroscopic redshift was used when available and the pho- 
tometric redshift was considered in the other cases (when- 
ever it was used our methods took into account its larger 
relative error) . Various photometric redshift catalogues have 
been constructed from SDSS data, and we compared these 
against spectroscopic redshift catalogues in different magni- 
tude ranges to choose which to use. We selected the pho- 
tometric redshift provided in the co l umn ph.otozcc2 of the 
schema Photoz2 (jOvaizu et alj|2008l : ICunha et aLll2009T ) for 
galaxies with a magnitude fainter than 17.77 because its dis- 
tribution with respect to the redshift was clearly the most 
accurate in the magnitude range 17.77 < r < 22.0. For 
galaxies without spectra and r < 17.77, we used the pho- 
tometric redshift provided in the column z of the schema 
Photoz, which seemed to agree better with spectroscopic 
redshifts for galaxies at these magnitudes. 



and e) passive galaxies. The types LINER, Seyfert and TO 
are aggregated and considered as the optical AGN. It is im- 
portant to notice that the catalogue does not include galax- 
ies that present broad emission lines as Seyfert 1. In practi- 
cally all of the galaxies (99.9 per cent) an [O ill] emission line 
(at 500.7-nm) with a luminosity of 10 6 ' 5 Lq would be de- 
tectable above the noise in our sample. To obtain a complete 
unbiased sample of active galaxies traced by optical emission 
lines, we used this cutoff for galaxies classified as optical 
AGN or as star forming (SF). The galaxies with [Om] lu- 
minosities below this threshold were included in the passive 
gala xies group. 

iBest fc Heckmanl (|2012h obtained the radio nuclear ac- 
tivity classification for galaxies in the MPA-JHU catalogue. 
The me thod followed the techniques presented in Best et al. 
(2005b) and used radio-continuum data from the National 
Radio Astronomy Observatory (NRAO) Very Larg e Ar- 
ray (VLA) Sky Survey (NVSS; ICondon et all [1998 ) and 
the Fain t Images of the Ra dio Sky at Twenty centimetres 
(FIRST; iBecker et al]|l995h databases. The separation be- 
tween "low-excitation" and "high-excitation" radio galaxies 
(LERG and HERG respectively) is also given. The radio- 
continuum luminosity limit at z = 0.1, corresponding to the 
flux density level (5 mjy) achieved in the cross-matching 
process, is Li.4GHz « 10 23 W m -2 Hz -1 . We consider as ra- 
dio AGN galaxies those that are classified as a radio AGN by 
IBest fc Heckmanl (|2012r ) and have a radio luminosity above 
this threshold. 

Table [1] shows a breakdown of the number of galaxies 
in each sample depending on their classification. 



3 ENVIRONMENTAL PARAMETERS 

We chose three different parameters to trace different as- 
pects of the environment and the interaction level. We de- 
fined a local density and a tidal forces estimator and ob- 
tained an estimation of group richness from the literature. 

3.1 Density estimator 

The density estimator (also called local density estimator) 
traces the density of companions around the target galaxy 
which is related to the dark matter halo density and is 
usually defined by the distance to the 5th or 10th nearest 
neighbour (e.g. lDresslerlll980l ; iMiller et al1l2003l ; lBestJl2004l ; 
IVerlev et al.ll2007r ). It can be defined as 

^^V^)=^4^)' 

where is the projected-distance (in Mpc) to the fc th near- 
est neighbour. It is possible to find different definitions of 
this parameter in the literature. I n some case s, a cylindrical 
volume or a surface density (e.g. Best 2004) is used. Note 
that the main factor that drives the parameter is the ratio 
fc/rfc and these different definitions produce only a scaled or 
shifted version of rjt. 

We dynamically obtained a list of potential compan- 
ions around each of our galaxies within a given radius. We 
considered as companions galaxies with a difference in spec- 
troscopic redshift lower than 0.01 (\z c — z t \ ^ 0.01) corre- 
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sponding to Av < 3000 km s 1 and with luminosities in r- 
band between 10 9 86 and 10 120 L . The lower limit of 10 9 ' 86 
Lq corresponds to the luminosity limit of the main spectro- 
scopic sample at redshift 0.1. Once selected, all the potential 
companions were supposed to be at the same distance as the 
target galaxy and the k- correct ions were comp uted using the 
analytical method of lChilingarian et alj (|2010T l. We used the 
distance to the 10th nearest neighbour. 

T he main spectrosco pic sample is ~94 per cent com- 
plete (|Strauss et al.l I2OO2T ) . The spatial distribution of the 
few galaxies with similar magnitudes that are not covered 
by the spectroscopic survey could affect the accuracy of the 
environmental estimation for some galaxies. However, given 
the size of the sample under study this is unlikely to be im- 
portant unless it introduces systematic biases. One indicator 
of a systematic bias in the local density is the presence of 
a dependence with redshift. If the number of companions is 
under- or over-estimated by our method for a given redshift 
range, a monotonic trend in the mean values of the esti- 
mator should be expected. In Fig. [1] the dependence of the 
local density (upper panel) with redshift is shown. There is 
no strong trend with redshift for the parameter. The varia- 
tion of the mean with redshift is less than 0.84 per cent with 
respect to the variation range of the local density. 



3.2 Tidal forces estimator 

The tidal estimator traces the relative tidal forces exerted 
by companions wit h respect to the in ternal binding forces of 
the target galaxy l|Verlev et al.ll2007r i. The definition is 

where D is the estimated diameter of the target galaxy and 
R is its radius, d is the distance between the target and the 
companion, Lr is the corrected luminosity in r-band and 
M is the mass of the galaxy. The luminosity in r-band is 
used, in this case, as a proxy for the mass of the galaxy. To 
broadly match the definition of diameter used in the litera- 
ture (projected major axis of a galaxy at the 25 mag/arcsec 2 
isophotal level or I Verlev et aill2007l l. we estimated the 
radius by using the Petrosian radius containing 90 per cent 
of the total flux (rgo; provided in the SDSS catalogue) scaled 
by a factor 1.46 (see Argudo-Fernandez et al. in prep.) 

For the tidal forces estimator, the inclusion of the pho- 
tometric galaxies (see Sect. 12. 3[) is critical, both because 
the SDSS fibre allocation process may not allow close com- 
panions to be targeted, and because the tidal forces are 
often dominated by nearby galaxies of lower luminosity. 
The first step before computing the tidal parameter was 
therefore to select which potential companions were actu- 
ally associated to our target galaxy rather than being fore- 
ground/background galaxies. For each potential companion 
galaxy without spectroscopy, we estimated the probability 
for it to be associated to the target galaxy by chance, and we 
used this probability to discard spurious companions. Specif- 
ically, for a hypothetical target at a redshift z, we computed 
the sky density of galaxies within the SDSS survey with an 
r-band magnitude brighter than a limit (rn m ) and with the 
difference between their photometric redshift and the red- 
shift of the target being less than N times the error of the 
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Figure 1. Dependence of the environmental parameters on the 
redshift. The relations of the local density (upper panel), tidal 
forces (middle panel) and cluster richness (lower panel) estimators 
with respect to the redshift are shown. The values for the target 
galaxies are shown as a grey-scale density cloud and dots for the 
outliers. Fifteen redshift bins where used. The mean (red solid 
line), its error (red dotted line), and the median (blue dashed 
line) for each redshift bin are plotted. The underlying black solid 
line marks the mean value for the whole sample. A random offset 
(ranging from —0.5 to 0.5) was introduced in the cluster richness 
parameter (n) to allow a proper visualization of the plot. 



photometric redshift. Hence, we obtained the surface density 
of galaxies as a function of z, ru m and N. For each poten- 
tial companion, the expected surface density of companions 
with its properties (r < r c , z — z t and N < z ^ zt ) was 
computed (the sub-indexes t and c denote target and com- 
panion respectively, z v is the photometric redshift and Az p 
is its error). If the expected number of companions within 
the area defined by the distance between the companion and 
the target galaxy is less than 0.05, we considered the can- 
didate galaxy as a genuine companion. Potential compan- 
ions with spectroscopy were selected if \z c — Zt\ ^ 0.01. The 
selected companions (photometric and spectroscopic) were 
presumed to be at the same distance as the target galaxy 
and k-corrections were applied. 

The search radius for tidal companions was limited to 
3 h _1 Mpc around the target (~ 4.29 Mpc using the cos- 
mology especified in Sect.[l}. We empirically checked that 
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increasing the search radius above this limit produced no ap- 
preciable change in the tidal parameter. If the relative force 
(tidal force over internal binding force) exerted by an indi- 
vidual companion is defined as qt,i = (Lrt/Ln)(2Rt/di i t) 3 
such that Q t = log(^~\ qt,i), the contribution of companions 
above the 3 h^ 1 Mpc limit was always q tA ^ 10 -6 , while the 
mean for closer companions is qt,i ~ 10 -2 . After the rejec- 
tion of spurious candidates, the tidal parameter was com- 
puted. In the rare case that no companions were selected 
(N = 684), an upper limit of Q = —6 was assigned to the 
galaxy. A galaxy ~ 100 times more massive than the target 
would be required to produce this effect at a distance of ~ 3 
h~ Mpc, our search radius limit. 

As a test for systematic biases, the middle panel of 
Fig. [T] shows the dependence of the tidal forces estimator 
with the redshift. In this case, the variation of the mean is 
less than 0.37 per cent of the variation range of the param- 
eter. 

It is important to note that the tidal estimator uses the 
projected two-dimensional distance as the distance between 
the target and the companion. This will lead to a general 
overestimation of the tidal parameter with respect to its ac- 
tual value. The overestimation can be more severe in regions 
with a high density of galaxies. However, as the projected 
distance depends only on the (random) line of sight with re- 
spect to the observer, we do not expect a systematic bias to 
be introduced by this effect. The expected effect is a possible 
decrease in the strength of the correlations between the tidal 
parameter and other properties, due to a degree of random 
scattering in the tidal paremeter measurement. 



3.3 Cluster richness estimator 

The last parameter is an estimation of the number of galax- 
ies in the cluster or gro up to which our target galaxy be- 
longs. iTago et al l (|2010h extracted groups of galaxies from 
the SDSS DR7 using a modified friends-of-friends (FoF) al- 
gorithm. A flux-limited and three volume-limited samples 
were provided. The groups of galaxies were detected and 
a value (n) was assigned to each galaxy corresponding to 
the number of galaxies in its own group, the cluster rich- 
ness. The number density of groups detected by FoF tech- 
niques depends strongly on redshift for a flux-limited sam- 
ple, hence, a volume-limited sample is required to avoid any 
biases. The volume-limited sample which covers all our red- 
shift range is the DR7-20. The cutoff is M r - 5 log(h) = -20 
(M r « -19.2) which covers from z = 0.02595 to z = 0.11012 
given the magnitude limits of the SDSS. Therefore, in our 
analyses, the reduced sample is composed of the galaxies that 
are bright enough in r-band and have an assigned value for 
the cluster richness. The value of the cluster richness ranges 
from 1 to 52. 

The dependence of cluster richness with redshift is 
shown in the lower panel of Fig. [T] The variation of the 
mean is less than 3.16 per cent of the variation range of 
the parameter. The presence of a few clusters with a high 
number of galaxies cause this slightly higher value with re- 
spect to the other parameters. The low values found for the 
variance of the mean of the three environmental parameters 
with respect to the redshift make us consider that they are 
not biased. 



Table 2. Principal Component Analysis results. 







Density 


Tidal 


log(n) 


VarP_ 


whole 


PCA1 


0.707 


0.707 




57% 


sample 


PCA2 


-0.707 


0.707 




43% 


reduced 


PCAlr 


0.661 


0.351 


0.663 


53% 


sample 


PCA2r 


-0.254 


0.937 


-0.242 


31% 




PCA3r 


-0.706 


-0.008 


0.708 


16% 



a Variance carried by the PCA component. 

3.4 Principal component analysis 

The three presented parameters provide three complemen- 
tary measures of the environment. The distribution of the 
local density and the tidal forces estimator is shown in the 
panel a of Fig. [2] Galaxies located at higher densities present 
medium to high tidal values but the highest tidal values are 
found for galaxies located in regions with medium to low 
densities. These galaxies are sometimes labeled as "field" 
galaxies according to their local density but they are clearly 
interacting. Hence, it is important to notice that the concept 
"field" and "isolated" are not equivalent. The panels b and 
c of Fig. [2]show the relation of the cluster richness estimator 
with the local density and the tidal parameter respectively. 
Each value of the cluster richness spans over a wide range of 
tidal and density values. In the case of the local density, it 
is clear that richer clusters present higher values of the local 
density. 

We performed a Principal Component Analysis (PCA) 
to break the parameters down into the independent physi- 
cal processes by removing the possible correlations between 
them. For example, part of the relation between the den- 
sity and the tidal parameter was expected to be produced 
by projection effects, with generally higher tidal values for 
galaxies that are located in denser environments. Some de- 
gree of relation between the group richness and the density 
was also expected, as found. 

We obtained two PCA components by combining the 
local density and the tidal parameters for the whole sam- 
ple and three components for the reduced sample where we 
added the group richness (the suffix "r" was added to distin- 
guish these principal components; we entered log(n) instead 
of n as the input of the PCA). The PCA components are 
shown in Table [2] and the normalization values applied to 
the environmental parameters are shown in Table The 
effect of the two-dimensional projection on the tidal param- 
eter was at least partially mitigated after considering the 
correction introduced by PCA2 with respect to the density 
(and cluster richness). 

As it is shown in Table [2] PCA1 and PCAlr are com- 
posed of a positive contribution of all the environmental pa- 
rameters. PCA2 results from a positive contribution of the 
tidal parameter and a negative contribution of the density 
estimator, with PCA2r also having a negative cluster rich- 
ness contribution. This similarity leads us to consider that 
PCA1 and PCAlr trace practically the same fundamental 
property, as do PCA2 and PCA2r. We later checked that 
the results obtained using either PCA1 and PCA2 or PCAlr 
and PCA2r were similar, hence, we will show the numbers 
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Figure 2. Environmental parameters. Panels a and d show the distribution for the whole sample and the rest of the panels for the 
reduced sample. The values for the target galaxies are shown as a grey-scale density cloud and dots for the outliers. A random offset 
(ranging from —0.5 to 0.5) was introduced in the cluster richness parameter (n) to allow a proper visualization of the plots. 



obtained using the tidal, density, PCA1 and PCA2 for the 
whole sample and the cluster richness and PCA3r for the re- 
duced sample. The relations between the PCA components 
are shown in Fig. [2] 

We interpret the PCA components as follows: 

(i) PCA1 and PCAlr follow the direction of the increase 
of the density (and the cluster richness) and also the tidal 
forces. We consider it to trace the overall interaction level 
and environmental nature of a galaxy. 

(ii) PCA2 and PCA2r are mainly driven by the differ- 
ence between the density (and cluster richness in the case 
of PCA2r) with respect to the tidal forces. A higher value 
traces higher one-on-one interactions while a lower value 
traces galaxies that are relatively isolated for its overall en- 
vironment. 

(iii) PCA3r is practically not affected by the tidal force 
and is driven by the difference between the density and the 
cluster richness. In this case, a galaxy within a rich group 
but with low local density values has high PCA3r values and 
vice-versa. Hence, it would probably trace scatter in den- 
sity versus cluster richness relation, being high where densi- 
ties are low for that particular n (e.g. cluster outskirts) and 
low in particularly locally-compact structure (e.g. compact 
groups). 

In Fig. O the relations between the PCA components and 
the environmental parameters are shown. PCA parameters 



Table 3. Principal Component Analysis normalization values. 





Parameter 


Mean 


Standard deviation 


whole 


Tidal 


-2.0275 


1.1952 


sample 


Density 


-0.5082 


0.7254 


reduced 


Tidal 


-1.8951 


1.1788 


sample 


Density 


-0.4865 


0.7100 




log(n) 


0.1708 


0.2748 



were computed applying the values presented in Table [2] to 
the normalised version of the environmental estimators ob- 
tained using the parameters presented in Table [3] The final 
values of the environmental estimators for all galaxies in the 
sample are presented in Table [4] 



4 RESULTS 

Our aim was to compute the relative fraction of each type 
of AGN with respect to the environmental and PCA pa- 
rameters. The strong dependence of the prevalence of AGN 
with the mass of the h ost galaxy in the case of o ptical 
jKauffmann et al.l2003bh and radio l|Best et alj2005al ) AGN 
has to be considered. Hence, to begin with, we divided the 
sample in four bins with masses between log(M) = 10 [Mq] 
and log(M) = 12 [M ] for radio AGN and six bins with 
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Table 4. Environmental parameters. The full version of the table can be found in the online version of the journal. 



plate 


mjd 


fiberlD 


Density 


Tidal 


log(n) 


PCA1 


PCA2 


PCAlr 


PCA2r 


PCA3r 


Mass 


Class.^ 


267 


51608 


1 


-1.318 


-2.881 




-1.294 


0.284 








9.62 


— S- 


286 


51999 


1 


-1.159 


-1.904 


0.000 


-0.562 


0.708 


-1.041 


0.384 


0.229 


10.63 


— AT 


291 


51928 


1 


-0.497 


-2.878 




-0.492 


-0.514 








10.25 




297 


51959 


1 


-0.943 


-2.546 




-0.730 


0.116 








10.32 




299 


51671 


1 


-0.189 


-2.570 


0.000 


-0.010 


-0.632 


-0.336 


-0.491 


-0.731 


10.77 


—AS 



a Classification code. The first pair of characters refer to radio activity and the second pair to optical activity. First 
character, radio AGN classification: R Radio AGN; - Not a radio AGN. Second character, type of radio AGN: H HERG; L 
LERG; - Non classified. Third character, optical AGN classification: A optical AGN; S star forming nuclei; - Passive. Fourth 
character, type of optical AGN activity: S Seyfert; L LINER; T Transition object; - Not an optical AGN. 
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Figure 3. Relation between the environmental parameters and 
the PCA. Panels a, b, d and e show the distribution for the whole 
sample and the rest of the panels for the reduced sample. The 
values for the target galaxies are shown as a grey-scale density 
cloud and dots for the outliers. A random offset (ranging from 
—0.5 to 0.5) was introduced in the cluster richness parameter (n) 
to allow a proper visualization of the plots. This random offset 
was also propagated to the PCA components that involve the 
reduced sample for the same reason. 

masses between log(M) = 10 [M©] and log(Af) = 11.8 [M©] 
for optical AGN; the AGN fractions as a function of envi- 
ronmental parameters are shown for each of these mass bins 
in Fig. |4] The monotonic increase of the fraction of active 
galaxies with respect to the mass of the host is clearly visi- 
ble. This increase is steeper for radio AGN than for optical 
AGN. For radio AGN the trends of AGN fraction with en- 
vironmental parameters are very similar for the four mass 
slices. In the case of optical AGN the trends are consistent 
for the density, PCA1 and PCA2 but for the tidal estima- 
tor, the cluster richness and PCA3r the results are not so 
clear. The bin of masses between log(M) = 11.5 [M©] and 
log(Af) — 11.8 [Mq] tends to not to agree with the others, 
probably because of the low number of galaxies within this 



mass range (389 galaxies in the reduced sample or 432 galax- 
ies in the whole sample which represent 0.4 or 0.2 per cent 
of the total number of galaxies considered respectively) . 

Although the overall trends could be inferred from 
Fig. |4l methods that allow an aggregation of the informa- 
tion of the different mass slices and hence a quantitative 
measurement of the effect of the environmental or PCA pa- 
rameter are preferable. With the aggregation, a better signal 
to noise ratio can be obtained and the quantification allows 
to draw more robust conclusions. Therefore, we first defined 
an aggregated ratio that joins the trends found in the differ- 
ent slices, and then, a statistical test that takes into account 
the effect of the mass (as a possible confounding factor) was 
applied. 

4.1 Stratified study 

To discard the effect of the mass on the fraction of AGN, we 
performed a stratified study of the ratio of AGN with respect 
to each environmental parameter. The sample was sliced in 
several mass strata (defined by the subindex i) and also in 
several bins of the environmental parameter considered (de- 
fined by the subindex j). This two-dimensional grid was 
fixed, and the following values were counted for each bin: a) 
riij; the total number of galaxies in a bin and, b) n a ,i,j', the 
total number of galaxies with a given type of nuclear clas- 
sification in a bin (we will use the term AGN from now on, 
although it can refer to any type of nuclear activity includ- 
ing passive galaxies) . The probability of a galaxy to harbour 
an AGN in each mass stratum is p a ,i = y] . n a ,ij/)] . riij. 
If there is no environmental dependence of the AGN frac- 
tion, then an estimated value of the number of AGN in each 
bin of the grid can be computed as n a _est,i,j = p a ,iTii t j. This 
estimate accounts for the dependence on mass, so any re- 
maining effect is due to environment Q Hence, a trend of the 
prevalence of AGN with respect to the environmental pa- 
rameter, that takes into account the effect of the mass, can 
be extracted by comparing n a _est,i,j and n a> i t j. The general 
trend with respect to the environmental parameter can be 
found by folding into the mass axis: n a _ est j = J^. n a _est,i,j, 

1 Notice that we remove largely the mass effects but small resid- 
ual trends can remain if e.g. AGN and all galaxies are differently 
distributed within the defined mass bins. 




Dens. Tidal log(n) PCA1 PCA2 PCA3r 



Figure 4. Fraction of optical or radio AGN galaxies, segregated in mass slices, with respect to the different environmental parameters. The bins with less than 20 galaxies (large error 
bars) were removed from the plots. Filled regions mark the ±f error zone. 
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n a ,j = ~^2,i n a,i,j and rij = X^ ni J- Then, the probabilities 
of harbouring an AGN with respect to the environmental pa- 
rameter are p a ,j = n a j/nj while the estimated value (taking 
account of the mass distribution) is p a _est,j = n a _ est ,j /rij . 
The errors on these values can be calculated following a 
binomial distribution. The final trend, defined by the rela- 
tive fraction between the observed and the estimated prob- 
abilities, is fi = p a ,i/Pa-est,i- A value of fi - Afi > 1 or 
fi + Afi < 1 means that, for this range of the environmen- 
tal parameter, there is a significant increase/decrease of the 
observed fraction of AGN with respect to the expected one, 
after correcting for the effect of the mass. 

In Fig. [5] /, is shown for the different types of nu- 
clear activity and environmental parameters. Seven bins 
were used in the mass axis and 10 bins in the environmental 
parameter axis defining a detailed trend. A coarse general 
trend using two bins is also shown. These two bins were 
separated at the median of each environmental or PCA pa- 
rameter except in the case of the cluster richness for which 
galaxies with n = 1 were grouped in the first bin and the 
rest in the second one. We also divided the cluster richness 
in two bins using the median after discarding galaxies with 
n — 1 and found that the trends were the same. 

The trends with respect to the density are clear and 
monotonic. There is an increase with density in the preva- 
lence of passive galaxies and radio AGN and a decrease for 
all the optical AGN types and SFN. The trends are simi- 
lar for the cluster richness but slightly less acute. The tidal 
parameter shows a non-linear trend for optically selected 
types with the lowest AGN fractions at intermediate values. 
This may be related to correlations with other parameters. 
For example, the galaxies which are located in high den- 
sity regions tend to have medium tidal values (see panel a 
in Fig. [2]), and the strong trends found with respect to the 
density might therefore affect the relative ratio for galaxies 
with mid tidal values. This interpretation is supported by 
the fact that there is a direct clear relation between PCA2 
and the presence of all types of activity; this also demon- 
strates the power of the PCA analysis in separating the dif- 
ferent physical processes. The trends with respect to PCA1 
are similar to the density ones but lowered, especially for 
galaxies which present star formation (TOs and SFN). It 
is noteworthy that while radio AGN and passive galaxies 
show the same general trend with respect to density, they 
show opposite trends with respect to the tidal and PCA2 
parameters. The relations with respect to PCA3r are usu- 
ally non-linear (see F ig. [51); these type of n on-linear trends 
were also observed bv lPimbblet et all (|2012T ) for the fraction 
of optical AGN with respect to the distance to the center of 
a cluster. This non-linearity weakens the trends found using 
only two bins, making it difficult to draw strong conclusions. 
Nevertheless, some results are apparent. A lower fraction of 
optical AGN is generally found for lower values of PCA3r 
(high density smaller systems, e.g. compact groups). There is 
also a possible decrease of the fraction of LINERS and TOs 
towards higher values of PCA3r (e.g. cluster outskirts). 

4.2 Cochran-Mantel-Haenzsel test 

The Cochran-Mantel-Ha enzsel test (CMH; ICochranl 1 1954 
iMantel fc Haenszelll 19591 ) can be used to study the strength 
of the association between two bi-valued variables (odds ra- 



tio) and its significance after taking into account the effect 
of the possible confounding factors. In our case, we corrected 
for the effect of the mass. The first variable was the "expo- 
sure" to the environmental parameter and the second vari- 
able was the effect of this exposure on the triggering of cer- 
tain type of nuclear activity. We discern the "exposure" to 
the environmental parameter by separating galaxies above 
and below the median value for a given environmental or 
PCA parameter. In the case of the cluster richness, more 
than a half of the galaxies have an assigned n of 1. Hence, 
the two selected groups were galaxies with n = 1 and galax- 
ies with n ^ 2. 

We tested if the prevalence of galaxies with certain type 
of activity is independent of the environment or interaction 
traced by the belonging to one of the two groups described 
before (galaxies with high or low values of the environmen- 
tal parameter or PCA). We obtained the CMH ratio and its 
significance. The CMH ratio is the relative probability be- 
tween the groups with high and low environmental param- 
eters. A ratio compatible with 1 means that the prevalence 
of nuclear activity is independent of the environmental or 
PCA parameter or that a dependence cannot be found with 
the data used. A significant ratio (defined by the p- value), 
higher or lower than 1, indicates that this environmental or 
PCA parameter is increasing or decreasing the chance for a 
galaxy to harbour this type of nuclear activity. 

The results of the CMH statistical test are shown in 
Table \S\ and in Fig. [5] they are shown in a way that al- 
lows a quick visual comparison between them. We added in 
this case a separation of radio AGN into HERG and LERG 
types. An increase in the local density produces a decrease 
in the prevalence of optical AGN (LINER, Seyferts, TOs) 
and SFN. The effect is the opposite for radio AGN and pas- 
sive galaxies. When we separate radio AGN in HERG and 
LERG, they present opposite trends with respect to the lo- 
cal density: LERGs show a clear increase with density while 
HERGs appear to show a decrease like optical AGN (though 
uncertainties are large due to the small sample size). The ef- 
fect of the tidal parameter is not clear for passive galaxies 
and TOs (as expected from the non-linear dependence in 
Fig. [5]). A slight decrease of the fraction of LINERs and 
Seyferts with increasing tidal parameter is found while an 
increase is found for SFN and radio AGN. This might be 
caused by the non-linear relation due to the effect of high 
density galaxies explained before. The trends with respect 
to the cluster richness are similar to the trends with den- 
sity. The effects that depend on PCA1 are very similar to 
the ones driven by density except for SFN. Probably the de- 
creasing prevalence of SFN with density is compensated by 
the increasing prevalence with the tidal estimator to make 
the PCA1 relation flat. In the case of increasing PCA2, there 
is an increase in the AGN fraction for Seyferts, TOs, SFN 
and radio AGN and a clear decrease for passive galaxies. The 
trends for HERGs and LERGs are compatible here. PCA2 
reveals much stronger positive trends for optical AGN than 
tidal does, due to removal of density selection effects, i.e. 
using PCA we show the importance of the one-on-one inter- 
actions. The radio AGN trend with PCA2 is much weaker 
than that with the tidal parameter, indicating that much 
of what is seen with respect to the tidal parameter is just 
a density effect (which gets stronger going from density to 
PCA1). The trends with PCA3r generally tend to be weaker 
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Figure 5. Relative fraction of galaxies of one type, corrected for the effect of mass, with respect to the different environmental parameters. 
The blue thin solid lines mark the detailed trend of / with respect to each environmental parameter. The thick solid lines defined by 
two points show the coarse general trend of /. The colour of the thick colour line indicates whether a fiat trend is compatible within the 
error (orange) or not (increase - green; decrease - red). 



and more variable for optical AGN types, as expected from 
the non-linearities shown in Fig. [5] However, radio AGN 
and SFN have clearly increasing fractions with an increas- 
ing PCA3r. 



5 DISCUSSION 

5.1 Environment and interactions 

Of the parameters studied, mass is confirmed to be the main 
driving facto r for the triggering of an AGN, especially for 
radio AGN l|Best et all l2005al ; iKauffmann et al.1 12004T ) as 
shown in Fig[4] However, both local density and interactions 
also play a significant role on the prevalence of nuclear activ- 
ity. Fig. [6] and Table [5] show the clear relations between the 
prevalence of the different types of nuclear activity with the 



environmental parameters. The relations obtained show the 
clearly different and even opposite trends found for the AGN 
fraction depending on the environmental or PCA parameter 
used and the type of AGN studied. These differences may 
have led to previous discrepancies found in the literature. 

The fraction of optical AGN decreases toward higher 
values of the local density, PCA1 and the cluster richness, 
with the prevalence of AGN in the densest local environ- 
ments being a factor 2-3 lower (at fixed mass) than that 
in sparser environments. The similarity between the trends 
with respect to local density and cluster richness reflects that 
both parameter are closely related, as shown in the panel c 
of Fig. [2] with galaxies with higher n located at higher densi- 
ties. In rich environments there are several mechanisms that 
can remove or change the properties of the cold gas needed to 
fuel the star formation or AGN from galaxies: the stripping 
of diffuse gas around galaxies or strangulation, ram-pressure 
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Table 5. Statistical study results from the CMH test. For each cell of the table the first row shows the CMH 
common odds ratio and its error and the second row shows the statistical significance or p-value (i.e., the 
probability of the trend occurring by chance) measured by the CMH test. The hypothesis tested was that 
the observed nuclear activity type is independent of the environmental or PCA parameter. The typeface of 
the p- value depends on its value: bold if p < 0.01; bold italics if 0.01 ^ p < 0.05; regular if 0.05 ^ p < 0.1; 
and italics if p ^ 0.1 





Density 


Tidal 


log(n) 


PCA1 


PCA2 


PCA3r 


Passive 


1.358 ±0.061 


0.950 ± 0.043 


1.248 ±0.077 


1.143 ±0.052 


0.745 ± 0.034 


0.989 ± 0.055 




0.0000 


0.0296 


0.0000 


0.0000 


0.0000 


0. 7025 


LINER 


0.77 ±0.11 


0.84 ±0.13 


0.87 ±0.14 


0.80 ±0.12 


1.08 ± 0.16 


1.18 ±0.19 




0.0012 


0.0353 


0.1248 


0.0053 


0.3375 


0.0642 


Seyfert 


0.801 ±0.054 


0.898 ±0.061 


0.783 ±0.070 


0.807 ±0.055 


1.118 ± 0.076 


0.914 ±0.076 




0.0000 


0.0029 


0.0000 


0.0000 


0.0018 


0.0406 


TO 


0.807 ± 0.076 


1.08 ± 0.10 


0.84 ±0.10 


0.906 ± 0.086 


1.35 ± 0.13 


0.97 ±0.11 




0.0000 


0.1568 


0.0097 


0.0510 


0.0000 


0.6727 


SFN 


0.608 ±0.051 


1.39 ± 0.11 


0.787 ±0.100 


0.999 ±0.082 


1.84 ±0.15 


1.16 ±0.13 




0.0000 


0.0000 


0.0005 


0.9783 


0.0000 


0.0085 


Optical 


0.797 ± 0.042 


0.939 ± 0.050 


0.808 ± 0.056 


0.831 ± 0.044 


1.181 ± 0.062 


0.964 ±0.061 


AGN 


0.0000 


0.0253 


0.0000 


0.0000 


0.0000 


0.2790 


Radio 


1.55 ±0.21 


2.77 ±0.52 


2.00 ±0.28 


2.13 ±0.36 


1.31 ± 0.18 


1.20 ±0.16 


AGN 


0.0000 


0.0000 


0.0000 


0.0000 


0.0003 


0.0113 


LERG 


1.63 ±0.23 


2.96 ± 0.59 


2.07 ±0.30 


2.21 ±0.38 


1.30 ±0.18 


1.19 ±0.16 




0.0000 


0.0000 


0.0000 


0.0000 


0.0004 


0.0185 


HERG 


0.55 ±0.27 


1.01 ± 0.52 


1.07 ±0.56 


1.09 ±0.57 


1.47 ±0.75 


1.81 ±0.97 




0.0720 


0.9881 


0.8496 


0.8141 


0.2897 


0.1248 
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Figure 6. Cochran-Mantel-Haenzsel ratios for the different types of nuclear activity with respect to the environmental parameters. The 
error bars mark the 95 per cent confidence interval. The shape of the symbol indicates the confidence level of the hypothesis tested (the 
probability of harbouring this type of activity is not affected by the environmental parameter) measured by the p- value: circle if p < 0.01; 
square if 0.01 p < 0.05; triangle if 0.05 ^ p < 0.1; and cross if p ^ 0.1. 



stripp ing and galaxy harassment (see Ivon der Linden et all 
l20ld . and references therein). The stripping of cold gas and 
its warming can explain the decrease of the prevalence of 
optical activity (AGN and star formation) found in denser 
and richer environments. 



On the other hand, one-on-one interactions, traced by 
PCA2, produce an enhancement on the prevalence of SFN, 
TOs, Seyfert and radio AGN and a decrease of the frac- 
tion of passive galaxies. This agrees with an scenario where 
galaxy interactions can fuel gas to the central region of 
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the galaxy 



lHaan et alj 



Shlosman ct al. Hl99d;lBarnes fc Hernquistlll99ll; 



2009; Liu et al]|201ll ). However, notice that this 



is not the only mechanism that triggers nu clear activity (e.g. 
ISilverman et alJl20lj ; ISabater et alj|2012h : the AGN preva- 
lence only increases by a factor ~ 2 between galaxies with 
no nearby neighbours to those undergoing the strongest in- 
teractions. For optically selected nuclear activity, the values 
of the CMH ratios with respect to PCA2 (see panel e in 
Fig. E]) decrease in this order: SFN, TOs, Seyfert, LINERs 
and passive galaxies. This relation can be compatible with 
the existence of a time delay between the tidal interaction 
and the setting of the different types of activity with an 
interaction inducing first star formation, then a Seyfert-like 
AG N, and then a LINER, before returning to a passive state 
fe.g. lLi et al]|200Sl ; lDarg et al.ll20ld : IWild et al.|[2010T l. 

The slightly decreasing trend of AGN prevalence found 
for optical AGN (both LINERs and Seyferts) with respect to 
the tidal parameter (see panel b in Fig. O may be produced 
by an overestimation of the tidal parameter in regions with 
high density of galaxies due to the use of projected distances 
(as dicussed earlier). Another explanation is that galaxies 
with a high tidal parameter located in dense environments 
lack the supply of gas required to fuel their nuclear activity 
despite their ongoing one-on-one interaction (see next Sec- 
tion). When this effect is corrected for using the correlation 
found with the PCA by the use of PCA2, the increase of the 
prevalence of all types of active galaxies with PCA2 becomes 
conspicuous. This suggests that PCA2 may trace better one- 
on-one interactions than the uncorrected tidal estimator, or 
at least, it may better trace the sort of one-on-one interac- 
tions that can trigger nuclear activity (i.e. those that fuel 
gas into the nucleus). 

Interactions and environmental density are complemen- 
tary factors that must be taken into account to get a com- 
plete view of the external phenomena that affect the trig- 
gering of an AGN. The trends found can be explained if the 
different parameters trace different physical processes that 
affect the triggering of nuclear activity. 



5.2 AGN types 

Optical AGN and HERG are thought to be fuelled in dif- 
ferent ways than the prevalen t population of radio AGN 
fLERG; lBest fc Heckmanll2012l ). The first type of AGN are 
fuelled at relatively high Eddington rates in radiatively ef- 
ficient accretion discs while the second type is fuelled at 
relatively low Eddington rates in radiatively inefficient ac- 
creting flows. While the first type requires a plentiful sup- 
ply of cold gas, the lower accretion-rate requirements of the 
LERGs could be achieved by accreting warm gas located in 
the haloes surrounding the galaxy, cluster or group. 

The local density parameter can trace the content and 
physical condition of the gas. The lack of cold gas in clus- 
ters of galaxies in the lo cal Universe is well known (e.g. 
iGiovanelli fe Havneslll985l ). and the gas that can be f ound 
at higher densities is warmer fsee lCattaneo et al.ll2009h . If a 
plentiful supply of cold gas is required to fuel optical AGN 
and SFN, then these types of activity will be suppressed 
in dense environments but enhanced by one-on-one inter- 
actions, as discussed above. On the other hand, the accre- 
tion of warm gas in would explain the higher prevalence of 
radio AGN of LERG type in denser environments. The re- 



sult that LERGs are enhanced by one-on-one interactions, 
however, argues against the simplest models in which all 
LERGs are fuelled from the warm gas haloes around mas - 
sive galaxies, groups or clusters fe.g. lHardcastle et alj|2007h . 
It suggest that some LERGs are fuelled by other mecha- 
nisms, and supports arguments that the distinction between 
radiatively efficient and radiatively inefficient AGN is more 
likely to be a function of Eddinton-scaled accre tion rate (e.g. 
iMerloni fc Hein j|200d iBest fc Heckmarj|2012l ) than the ori- 
gin of the fuelling gas. In this picture, most galaxy interac- 
tions provide gas to the nucleus at a sufficiently high rate 
that a radiatively efficient (optical) AGN is produced, but 
some interactions (or at some times during an interaction) 
the gas flow to the nucleus is slower and results in a LERG. 
These LERGs may be distinguishable fro m the cooling- 
fuelle d LERGs in terms of their colours (e.g. Ijanssen et al.l 
l2012j ). 

The low numbers of HERGs make difficult to draw sta- 
tistically significant results from them, but their prevalence 
is statistically higher at lower than at higher local densities 
(see panel a in Fig. [6]). This is opposite to the trend found 
for LERGs but the same as that for optical AGN. Indeed, 
all the environmental trends for HERGs are more similar to 
optical AGN than to the rest of radio AGN (LERGs). This 
suggests that the physical processes that trigger HERGs are 
similar to that of optical AGN. Despite the similar pow- 
ering mechanisms and environment trends, it is important 
to note that HERGs do show significant differences with 
respect to Seyferts, most particularly that they are prefer- 
entially found at higher stellar masses (median values for 
the whole sample: 10 10 ' 88 vs. 10 10 ' 56 M ; reduced sample: 
10 10 - 98 vs. lO 10 ' 74 Ms). 

The mechanis ms that power LINE R galaxies are a mat- 
ter of controversy l|Kewlev et al.ll20p6|V Some LINERs are 
certa inly powered by accretion fe.g. iGonzalez-Martm et al.l 
2006), but there are strong arguments that a fraction of 
them could be powered by low-mass evolved stars (e.g. 
ICid Fernandes et al.ll201lh . In our analysis we found that 
the LINERs show very similar trends to those of the rest of 
the optical AGN population (although some differences can 
be observed in PCA3r). This may be caused by the limiting 
luminosity that we apply to Lront] > which will select only the 
most lum inous LINERs which are m ore likely to be powered 
by AGN (|Cid Fernandes et al.ll201ll); these may behave di f- 
ferently from low-luminosity AGN (|Schawinski et al] l2010). 
This may cause the trends of our LINERs to be similar in 
general to those of Seyfert galaxies. 



5.3 Further considerations 



In ISabater et all ||2012D . no evidence of a difference 
was found between the prevalence of AGN in iso- 
lated galaxies (AMIGA sam ple of isolated galaxies; 
IVerdes-Montenegro et al1l2005f) or in a sample of compact 
groups ( Martinez et alj 2010t ) at a fixed stellar mass and 
morphology. This result can be explained if the increment 
expected in compact groups due to higher one-on-one inter- 
actions is broadly compensated by their higher local den- 
sity and deficiency of cold gas (|Verdes-Montenegro et ahl 
2001). However, an interpretation of t he differences found 
for star formation in the literature (e.g. ISulentic et aTll200ll : 



14 J. Sabater, P. N. Best and M. Argudo- Fernandez 



iDurbala etTafl 120081 ; iMartinez-Badenes et ail 120121 ) is not 
straightforward. 

There are additional considerations that could be taken 
into account to further advance this study. Firstly, mor- 
phology could play an important role on nuclear activity 
|Hwang et al J 12013 ). and its effect as a confounding factor 
could also be estimated. Secondly, lower emission line lu- 
minosity systems could be probed to investigate the more 
typical LINER population. To allow the discrimination be- 
tween genuine AGN-powered LINERs and retired galax- 
ies whose emission resembles that of weak AGN, alterna- 
tive diagnostic methods that take into account the differ- 
ences between these types of weak AGN would have to 
be considered ([Buttiglione et al.ll2010l ; ICid Fernandes et al.l 
l201ll ; ISabater et al.ll2012f) . Thirdly, it could be investigated 
whether different sub-populations of LERGs (e.g. by colour 
or luminosity) show the same environmental and interaction 
dependencies. This study will be extended in the future to 
clarify those points. 



6 SUMMARY AND CONCLUSIONS 

We presented a study of the prevalence of optical and radio 
nuclear activity with respect to the environment and inter- 
actions in a sample of SDSS galaxies. We aimed to determine 
the effect of these external factors on the triggering of nu- 
clear activity but taking into account the effect of possible 
confounding factors like the galaxy mass. 

The study is based on a sample of ~270000 galax- 
ies with a redshift between 0.03 and 0.1 drawn from the 
main spectroscopic sample of the Sloan Digital Sky Sur- 
vey (SDSS). Additional data like the stellar mass and the 
optical AGN classification were obtained from the MPA- 
JHU catalogue. The rad io AGN classification is described in 
iBest fc Heckmanl lj2012h . We defined a local density param- 
eter and a tidal forces estimat or and used a clu ster richness 
estimator from the literature (Tago et al. 2010) to trace dif- 
ferent aspects of environment and interaction. A principal 
component analysis was applied to the environmental pa- 
rameters to consider and remove the possible correlations 
between them. A stratified statistical study, which removed 
the effect of the galaxy mass, was applied in order to find the 
unbiased relation of the prevalence of the different types of 
nuclear activity with the different PCA and environmental 
parameters. The relations were quantified and their statisti- 
cal significance obtained using the Cochran-Mantel-Haenzsel 
test. 

From this study we found that 

(i) There is a strong dependence of the prevalence of nu- 
clear activity with density and interactions after taking into 
account the effect of mass. However, it is fundamental to 
distinguish and consider both local environment and one- 
on-one interactions because they can produce different and 
even opposite effects. 

(ii) We found a decrease of optical nuclear activity, in- 
cluding star formation, towards denser environments. This 
effect is probably driven by the lack of a cold gas supply at 
higher densities. 

(iii) On the other hand, an increase of the prevalence of 
radio nuclear activity is found towards higher densities. This 



can be explained by the warmer gas present at higher den- 
sities that is accreted at low rates in a radiatively ineffi- 
cient manner, triggering typical low-luminosity radio AGN 
or LERGs. 

(iv) There is an increase of the prevalence of SFN, optical 
and radio AGN with one-on-one interactions. The depen- 
dence for the different types of optical activity with respect 
to PCA2 traced by the CMH ratio follows this order: SFN 
(1.84), TOs (1.35), Seyferts (1.12), LINERs (1.08) and, fi- 
nally, passive galaxies (0.75); this is in accordance with the 
time sequence that has been argued for the different activity 
types. 

(v) The trend for HERGs with local density is statisti- 
cally different to that of the rest of radio AGN (LERG) and 
similar to that of optical AGN. This adds evidence to the 
link of HERG with high excitation AGN, like optical and 
X-ray AGN, that needs a supply of cold gas to be powered. 

Our results agree with a scenario where the type of AGN 
triggered is related to the presence and physical condition of 
the gas supply of the AGN, which in turn depends on both 
large-scale environment and local interactions. High exci- 
tation AGN (optical AGN and HERGs) and SFN require 
the presence of a plentiful cold gas supply to the nuclear 
regions. This gas is less likely to be present in denser envi- 
ronments where galaxies tend to be gas poor due to stripping 
or harassment. The prevalence of AGN is enhanced in the 
presence of one-on-one interactions which can funnel gas to 
the nuclear regions. These two effects may work in compe- 
tition against each other (one-on-one interactions are more 
likely in denser environments) and hence, the consideration 
of different aspects of the environment is fundamental to 
obtain a complete understanding of the mechanisms that 
may trigger the different types of AGN. In contrast, low- 
luminosity radio AGN of LERG type require relatively low 
gas accretion rates. This can be supplied by the cooling of 
gas from the hot haloes of groups and clusters accounting for 
the strong increase in the prevalence of LERGs with local 
density. However, the presence of a residual trend with one- 
on-one interactions after density effects are removed (PCA2) 
suggest that warm gas cooling is not the only mechanism by 
which LERGs can be triggered. 
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